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Abstract 



We consider the Navier-Stokes equation on a two dimensional torus with a random force, 
white noise in time and analytic in space, for arbitrary Reynolds number R. We prove 
probabilistic estimates for the long time behaviour of the solutions that imply bounds for 
the dissipation scale and energy spectrum as R — * oo. 

^ ' 

-i— > 

a 

1 Introduction 

X 

In two dimensions global existence and uniqueness of solutions of the Navier-Stokes equation is 
known for a large class of initial conditions and forcing, deterministic and random. In particular, 
for a bounded domain the unforced system has a finite dimensional attractor Q and this 
persists for a bounded finite dimensional force. 

In this paper we consider the Navier-Stokes equation with a random force, white noise in 
time and large scale in space and prove probabilistic estimates for the long time behaviour of 
the solutions. Our analysis is inspired by the recent paper by Mattingly and Sinai || who gave 
a conceptually simple proof of analyticity of the solutions of the 2d Navier-Stokes equation. 
We extend their analysis to the random case. 

We consider the stochastic Navier-Stokes equation for the velocity field U(t, x) defined on 
the torus T L = (R/2vrLZ) 2 : 

dU+((U- V)V -uV 2 V + Vp)dt = dF (1) 

where F(t, x) is a Wiener process with covariance 

EF a (s, x)Fp(t, y) = min{ S , t}C a(S (^) (2) 
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and C a/ 3 is a smooth function defined on the unit torus and satisfying d a C al 3 = 0. (0) rep- 
resents large scale forcing, the scale being the size of the box. ([I]) is supplemented with the 
incompressibility condition V ■ U = = V • F and we will also assume the vanishing averages 
over the torus: f T U(0, x) = = J T F(t, x) which imply that f T U(t, x) = for all times t. 
(0) implies the transport equation for the vorticity Q = d\U2 — dzU\\ 

dVL + ((U • V)fi - vV 2 n)dt = dG, (3) 

where G = d\F 2 — d%F\ has the covariance 

EG(t,x)G{s,y) = L- 2 mm{s,t}T(^) 

with T = -AtrC. 

It is convenient to change to dimensionless variables s.t. v and L become one. This is 
achieved by setting 

U(*,x) = £u(£t, ~x) , f2(t,x) = ^(jpt, ~x). 

Then u and u live on the unit torus and satisfy (HD and @ with i/ and L replaced by 1, and C 
and T replaced by 

Going to the Fourier transform u^it) = {2ir)~ 2 J T e ik ' x u;(i, x)dx with k e Z 2 we may write the 
enstrophy equation as 

dcj k = (-k 2 cj k + £ (k x + d/ k (4) 

iez 2 \{o,k} 

where k x 1 = kih — hk^ and {/k} are Brownian motions with /k = /-k an d 

Efk(s)fi(t) = min{s,t}5 k -i7k 

and we have used the relation u k = i ( ~ h ^ kl) oj-^. 

The dimensionless control parameter is the u injection rate, 

R = \ E 7k = | 7 (o) = |^r(o), 

kez 2 

that plays the role of Reynolds number in our model. We will be interested in the turbulent 
region R — > oo. We make the following assumption on the noise covariance: 

7k < CRe-^. (5) 

The coefficient of |k| is arbitrary, but we require exponential decay. The physically relevant 
case is the one with 7k 7^ only for a finite number of k with |k| of the order of unity. 
To state our main result, define the enstrophy 



$ = ^EKI 2 (6) 



k 

and fix numbers r > 1, a > 1 + r. Consider, for positive D, the norm 



\u\\ D = sup|w k ||kre^ Q|k| . (7) 

k 
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D will vary below, but r and a are fixed. The factor |k| r is useful technically (and was already 
used in 0). 

Theorem. Let \\w(0)\\d o < Dq < oo and $(0) = K < oo. Then, there exists a random 
function D t , D t < oo for all t, such that with probability 1, \\w{f)\\D t < Df . For any t > 
C(log D + log K) , and for D 2 > CR log R, 

Piob{\\u{t)\\ D <D a k $(t) < D 2 } > 1 - Ce' c ^. (8) 

Remark. Here and below, C (and c) are sufficiently large (small) constants, which may vary 
from place to place but that are uniformly bounded as R — > oo. The theorem says that with 
probability one ui(t, x) is analytic for all times, the dissipation scale is (up to a logarithm) 
> R~^ a and the energy spectrum 

e(Jfe) = A;" 1 jT i rfk^l^J 2 < CR Z k~ {2r+1) 

with = |k|, where r can be taken arbitrary close to 1 and a arbitrary close to 1 + r. These 
bounds hold for any fixed time and also for the average of these quantities over any fixed time 
interval. For example, using Jensen's and Chebyshev's inequalities, one derives from (§) 

Prob{i j* +T '\tu k (s)\ 2 ds > D 2a k- 2r e~ 2D ' ak } < Ce' c ^. 

Let us close this section with two comments. The first concerns the relationship of our 
model to the standard 2d turbulence picture || One considers ([!]) in infinite volume with 
the forcing as we do at spatial scale L, but not periodic, rather, for instance, having a smooth 
Fourier transform with compact support around L^ 1 . Then it is expected that a stationary 
state for Q emerges for which the energy spectrum e(k) = k~ x f sl dh f dxe th]l ' x EQ(x.)Q(0) has 
two scaling regimes 

6(A;)0C U-t k «L-i ) (9) 

refered to as the direct (enstrophy) cascade regime and the inverse (energy) cascade regime 
respectively. The scale rj is the "viscous scale" beyond which the e(k) decays more rapidly and 
it scales like . In particular, the total energy density f °° e{k)dk is infinite in the stationary 
state. This means that starting with say vanishing u at time zero, the energy density increases 
linearly with time and for the ensuing stationary state only the vorticity remains a well defined 
random field. One can also work in finite volume like in this paper by forcing the system in 
an intermediate scale r] « £ « L, provided the energy is absorbed by friction acting on 
the |k| ~ L^ 1 regime. This indeed is what one does in experimental |J and numerical || 
approaches. 

In our case the absence of the friction forces the energy to dissipate in the short scales too 
and the spectrum should be different from ([]). Our bound above is certainly far from realistic, 
but one would expect the e{k) to diverge as R — > oo. It would be very interesting to get hold of 
the direct and inverse cascade regimes, but certainly much more sophisticated ideas are needed 
than what are used in the present paper. 

The second comment concerns the uniqueness of the stationary state (the existence is stan- 
dard and follows from compactness and Lemma 1 below). In the case of Gaussian noise like as 
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we have there are two kinds of results in the literature regarding uniqueness. In |J one proves 
uniqueness, provided the noise is taken big enough in the ultraviolet, i.e. the 7k are taken 
to have a lower bound k~ a for a sufficiently small. This assures that ergodicity results from 
the action of the noise. However, such a noise is not what one is interested in the turbulence 
problem. The second result is for a smooth noise but viscosity large enough, i.e. in the 
nonturbulent regime. Then the Laplacean is the dominant term in equation (fl|) and the past 
is forgotten exponentially fast due to the viscous damping. 

In the turbulent regime of large R, the number N of modes that are not explicitely 
damped by viscosity goes to infinity as R — ► oo (we get an upper bound CR a for N). Never- 
theless, in the absence of noise, the enstrophy and thus u(t) tends to zero and this dissipativity 
should lead to uniqueness of the stationary state provided the noise is nonvanishing for these 
N modes (in the case of bounded noise, kicked at discrete times, uniqueness has been recently 
proven by Kuksin and Shirikyan UTOf). 

From the physical point of view, the rate of convergence to the stationary state that could be 
obtained solely due to the effect of the noise would not be realistic. As R — > oo, the relaxation 
time due to this mechanism would presumably grow superexponentially in R while in actual 
fact relaxation to stationarity should be due to the nonlinearity and should be much faster. 



2 Transition probabilities 

Define the region 

Uo = {to | < D a and $ < D 2 } (10) 

Then the basic proposition is 

Proposition. Suppose uj(0) E Ud- Then there are positive constants A and a, independent of 
R, such that 

Prob{cj(t) G U^— tD , Vt, < t < 1} > 1 - AR 2a e~ aI £ (11) 

It has a rather immediate 
Corollary. Suppose uj(0) G Ud and D' > D. Then 



V2 



Prob{a;(t) £ U^=- tD , , for some t G [0, 1]} < AR 2a e~ a ~ (12) 

Proof. Note that for D < D', Ud C Ud 1 - Thus u>(0) E Ur>>. Now the Proposition implies the 

claim. □ 

Proof of the Theorem. Consider the Markov chain with transition probabilities 

p(u, U) = Prob{tu(l) G U | uj(0) = uj}. (13) 
Let U n = Ud u where D 2 = 2a~ l R(\e) n and define 

Pm,n = SUp p(u, U°). (14) 



Since, by definition, y / 2e~ 1 D m < D n , for m < n + 1, the Corollary implies 

Pm,n < A'e~^ er = A'lT n (15) 



4 



for m < n + 1 and D 2 n > CRlogR (so that R 2a e~ a ^ < C, and we can take A' = C'A). 

By assumption, u(0) G Un for any N < oo such that D in the theorem is less than D N . 
Let p n (t) = Prob{cj(t) G Then 

Pn(^ + 1) < Prob{^(t) G t4+i}p„ + i in +p n +i(<) < Pn+l,n + Pn+1 

(*)■ (16) 

Suppose, inductively in t G N, that 

Pn(t) < Bn n (17) 
for n > N - t. Then, for n > N - 1 - 1, (TfJ) and © yield 

Pn(* + 1) < + 57T n+ i = Bn n 

provided we take B = A' (I — e~^ e ^) _1 (for t — 0, (|T7|) holds for any £> > 0). This completes 
the induction and shows that, with probability one, u(t) G U n , for some n, for all integer times. 
Moreover, since ( |T7D holds for all n when t > N = C (log Do + logK), this finishes the proof of 
the Theorem for integer times. The remaining times follow from the Corollary. □. 

3 Enstrophy bounds 

We prove a probabilistic analogue of the enstrophy balance: 
Lemma 1. Given $(0), for any t G [0, 1], 

Prob{$(i) > D 2 } < c e -*^ 2 -*(o)) 

Proof. Let x{t) = 2A(t)$(t) = A(t) Ek K| 2 - Then by Ito's formula (recall that £ k 7k = 2R 
and thus that 7^ < 2R, Vk): 

±E[e*] = £[(AA- 1 o:-2A]>>Vkr + A£7k + 2A 2 £7kKlV] 

k k k 

< ^[((AA- 1 - 2 + 4XR)x + 2XR)e x ] 

where E denotes the expectation taken over the /k's. We used the Navier-Stokes equation (H), 
|k| > 1, and the fact that that the nonlinear term does not contribute. Take now A(t) = gfee^ 1 ) 
so that AA" 1 = 1, AA" 1 - 2 + AXR < -\ and 2Ai? < \. So, 

±E[e x ] < E[(- - -x)e x ] <-- -E[e x ] 
m 1 1 — LV 4 2 2 4 

where the last inequality follows by using (1 — 2x)e x < 2 — e x . Thus, Gronwall's inequality 
implies that: 

E[e°®] < e-3e x(0) + 2 < 3e x(0) 

i.e. 

£[exp(|$(t)e')]<3exp(|$(0)), 

with c = ^-r- which yields the claim by Chebycheff 's inequality. □ 

This implies immediately the 
Corollary. Let D(t) = e~**D with D 2 = $(0), and let t h . . . ,t N G [0, 1]. Then, 

Prob{$0g < ^Ditn) 2 , Vn = 1, . . . , N} > 1 - CNe~ c ^. (18) 



4 Proof of the Proposition 

As usual, the stochastic equation (Q) is defined by the integral equation, 

u w (t) = e-* k2 ^ k (0) + fdse^ 2 £ (k x 1)|1|- + z k (t). (19) 

iez 2 \{o,k} 

where Zk is an Ornstein-Uhlenbeck process i.e. Gaussian with mean zero and covariance 

Ezjfyzfa) = 5 k ,-ii(e- ( *- s)k2 - e-^ k2 ) 7k 

Our strategy to prove the proposition is the following. We fix a short timestep r depending 
on D. By the Corollary of the previous section the enstrophy can be assumed to satisfy the 
required bounds at discrete times t n = nr. On the interval [0, r] we prove an existence and 
uniqueness result for (|T9|) in Lemma 3 by imposing a suitable condition on the smallness of 
the noise term z. At this point, the bound for ||w(t)||D(t) will not improve as claimed in the 
Proposition. However the enstrophy stays bounded and this information allows (Lemma 4) to 
improve the ||cj(t) Horn-bound. Repeating lemmas 3 and 4 on intervals [t n , t n+ i] the Proposition 
follows. 
Let 

T = 5D- 4a 



where 5 will be chosen below (see after (|29f) ), independently on D. We need the following 
standard result on the Ornstein-Uhlenbeck process: 

Lemma 2. Vk G Z 2 , Prob{sup t6[0jT] \z k (t)\ > Bt^} < Ce'*^^ 
This has the following simple consequence. Let Ad be the event 

{z \ Vk G Z 2 , sup \z k {t)\ < r*L>e-^}, (20) 
te[o,r] 

then Lemma 2 implies 

Prob^ D > 1 -Ce~ c ^-. (21) 



We now prove two lemmas. The first one, as we explained above, shows that the solution 
exists and that the solution satisfies the bounds of the Proposition over a short time interval. 
For this, let Yd be the Banach space equiped with the norm || • \ \d and 

X D = {coE C°{[0,t],Y d ) I |M| = sup \\to(t)\\D(t) < oo} (22) 

te[o,r] 

where 

D(t) = e'^D. 

Then we have, 

Lemma 3. Let z G Ad and suppose that ||o;(0)||d < D a and that $(0) < |-D 2 . Then the 
solution exists in Xd and moreover, 



(t)\\v2D(t) < (V2D(t)r , $(t)<2D(tf 



fortE[0,r]. 



The second lemma improves on these bounds: 

Lemma 4. Let z G Ad and suppose that \\u(0)\\ D < D a and that < 2D(t) 2 fort G [0,r]. 
Then\\tu(r)\\ D{T) < D(r) a . 

Proof of the Proposition. Let t n = nr. By the Corollary in the previous section fll8|), we 
may assume that $(i n ) < |D(t n ) 2 , for all n = 1, . . . , N, where iV = <5 _1 .D 4q - 1 with probability 

l-CD 4a e~ cJ £. (23) 

We can thus repeat Lemmas 3 and 4 on intervals [t n ,t n+ i], each time with probability (pT|). 
Hence, with probability bounded from below by (p3|) we deduce that ||^(i)II^D(t) 
and $(t) < 2D{t) 2 for all i G [0, 1], i.e. we have u(t) G U y ^—t D as required. By changing c and 
C, we can bound -j^- in ( |2"B| ) by the exponential, call a and A the new constants and obtain 
the claim of the Proposition. □ 

Proof of Lemma 3. Write equation ( |j"9"|) as 

uj = F(u) (24) 

where 

F k (v) = u;£(t) + fdse^ 2 ]T (k x = ^(t) + iV k (^)(t) (25) 

iez 2 \{o,k} 

and a>°(£) equals: 

u;°(t)EEe-* k2 ^ k (0) + z k (t). (26) 
Using ( p0[) and Zk(0) = 0, which imply (trivially) that \\z\ \ < D a , and 

e -tk" e -D-«|k| < e -D(t)-|k| (2?) 

which holds for t G [0, r], we have 

||cj°|| < 2D a . (28) 

We prove now that F is a contraction in the ball 

B = {v G X D : \\v-u°\\ < 1}, (29) 

provided the 5 in r = <5.D~ 4a is taken small enough (independently of D). To show that F 
maps B into itself, let i> G B. Then ||t> || < 2D a + 1 i.e. 

k(t)| < (2.0° + l)e- D(t) " Q|k| |kr r . (30) 

We must prove that 

\F k (v) -wj>(t)| = |JV k (t;)(t)| < e-^-Mlkr, (31) 



Vk G Z 2 and Vt G [0,r]. Inserting (^) and |k x 1| |1|~ 2 < |k||l| _1 in the second term of 
we get: 

|Mc(v)(t)| < (2 J D" + 1) 2 fdseW* Y, e- B W" l k - 1 le- D W" a l 1 l|k-l|- r |l|- r - 1 |k| (32) 

iez 2 \{o,k} 
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Then, using the bound 

£ |k-l|- r |ir r - x < C|k|- r , (33) 
iez2\{o,k} 

(since r > 1), the triangle inequality — |k — 1| — |1|<— |k| and 

|(s - t)k 2 < (e? as - e? at )\k\D~ a = (D(s)- a - D(t)- a )\k\, (34) 

which holds for < s < t < 1 and D large enough, one gets that 

\NM(t)\ < (2D a + l) 2 \k\C\k\- r e- D ^ a M C dse^ s - t)v? 

Jo 

= (2D a + l) 2 C|k|- r e- D W"'" |k| 2|kr 1 (l - e'^). (35) 

Since |k| _1 (l — e _7 * k ) < t^ < 5^ D~ 2a (|3l"l ) follows for 5 small enough (but independent of D). 
The contractive property is proven similarity. 

Combining the fact that the solution is contained in the ball (|29"D and the inequality 2D a + 
1 < (V2e~2D) a = (\/2D(t)) a (which holds, since a > 2, for t e [0, r] and D large enough) we 
obtain, 

< \\uj(t)\\D(t) < (V2D(t)) a (36) 

for t £ [0, r] and D large enough. 

To conclude we need to prove that $(t) < 2D 2 . By (|4|, |^), 

= < |(||^°(t)|| 2 + ||iVH(t)|| 2 ) 2 (37) 

By (0) and z k (0) = 0, the L 2 -norm of z(t) is bo unded by Cr^ D = Cd^D 1 - 2 *, the L 2 -norm 
of the first term in @ is bounded by | |w(0)|| < y/2$(0) < V3D and, using (pi]), the L 2 -norm 

of N(u>)(t) is bounded by (Skez 2 \{o} |k|~ 2r )^ = C (since r > 1). Thus, we obtain the claim 
provided D is large enough. □ 

Proof of Lemma 4. We note first that $(r) < 2D{t) 2 implies 

Mr) | < V2D(t) < ,D(r) Q e" D(TrQ|k| |kr r (38) 

provided |k| < D@, a > 1 + r/3 and .D is large enough. Hence, we only need to consider 
|k| > D 13 . Below, we take as (3 any number strictly larger than 1. 

We can now conclude the proof of the Lemma by using the following bound on the nonlinear 
term of the Navier-Stokes equation, which improves fl3"5| ) for k large enough: 

Lemma 5. Vk such that |k| > D 13 , andWt £ [0, t\, 

\N k (u)(t)\ < c(l - e -* tk2 ) J D a e _D( * ra|k| |kr r (39) 
where c can be taken small if D is large enough. 

Returning to the proof of Lemma 4, we have to prove the following bound: 

Mr) | < e~ iaT D a e~ D( - T) ~ a W\k\- r . (40) 
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We shall use ( [2"5| , |2"6j ) and bound each term. Using ||o;(0)||_d < D a and 

e -irk» e -U-[k| < e -%)"«|k[ (41) 



which is similar to (p7), we get, 

|e- rk2 cu k (0)| < e -* Tk? £) a e - IJ W- a l k l|k|-''. (42) 



-M < e -D(r)-|k|| k |-r (43) 



For z k (r), use ©, z k (0) = 0, and 



e 

for |k| large, to get 

W{t)\ < T?De-^e- D ^- a W\\L\- r . (44) 
Finally, we use (||) to bound N k (u)(r). Combining (g2|), (f44|), ©, we obtain g^) using 

e"i Tk2 + c(l - e-* Tk2 ) + r^D^e-^ < e-? ar . (45) 
Since r = 5D _4a , this last estimate holds for c small, |k| > and D large enough. 

Proof of Lemma 5. Consider first the case D* 3 < |k| < AD a , where A is a large enough 
constant (chosen below). We bound |k x 1||1|~ 2 < |k||l| _1 and split the sum in fl25"D into 

( E + E )K_i( s )ih(,)||k||ir 1 = s 1 + s 2 . (46) 

o^\i\<M Mk,|i|>M 
In the first sum, we bound, using Lemma 3, 

K_i(s)| < CD Q |k - l|- r < CD a |k|- r 
since |k — 1| > y|k|. Also, from Lemma 3 

(s)|| 2 = $(s)*< v^D(s) (47) 



so Schwartz' inequality yields 



E kOOllil-^vSiw E |ir 2 <^(iog|k|)- (48) 



o^|i|<H o#|i|<l£l 



Combining these two bounds we get 

£i < CD\k\{\og \k\)? D a \k\~ r . (49) 



For the second sum, we use |^>i(s)| < CD a |l| r , together with pTj ) and Schwartz' inequality 
to bound it by 

S 2 < CD\k\D a ( Yl |ir 2(r+1) )* < CD\k\D a \k\- r . (50) 

¥k,|l|>^ 

Inserting (f49|) and (0) to iV k (u;)(i) and performing the integral over time we get the bound 

\N k (u)(t)\ < (7D|k|- 1 (log|k|)*(l- e -* k2 )I} a |k|- r 

< Ce CA D^ (log D)^(l - e-^ 2 )D a e- D ^ a ^\k\- r (51) 



where we used D 13 < |k| < AD a and 

i < e -m-^\ e cA 

which holds since |k| < AD a . The claim of the Lemma follows, for D@ < |k| < AD a , since D 
is assumed to be large enough and we choose (3 > 1. 

Consider now the case |k| > AD a . Using the bound (|35|) , we get 

\N k (u)(t)\ < C\k\-\1 -e-^ 2 )D 2a e- D{t ^\k\- r < c(l - e^ t]l2 )D a e~ D ^ a ^\k\- r (52) 

by choosing A large enough (thus, we first choose A large so that fl52|) holds with c small enough 
for (|45|) to be true and then we choose D large so that the RHS of (|5lD is bounded by the RHS 
of (|39| ) with c small enough). □ 
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